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Abstract. The nature of dark matter (DM) particles and the mechanism that 
provides their measured relic abundance are currently unknown. In this paper we 
investigate inert scalar and vector like fermion doublet DM candidates with a charge 
asymmetry in the dark sector, which is generated by the same mechanism that provides 
the baryon asymmetry, namely baryogenesis-via-leptogenesis induced by decays of 
scalar triplets. At the same time the model gives rise to neutrino masses in the 
ballpark of oscillation experiments via type II seesaw. We discuss possible sources of 
depletion of asymmetry in the DM and visible sectors and solve the relevant Boltzmann 
equations for quasi-equilibrium decay of triplet scalars. A Monte-Carlo-Markov-Chain 
analysis is performed for the whole parameter space. The survival of the asymmetry 
in the dark sector leads to inelastic scattering off nuclei. We then apply bayesian 
statistic to infer the model parameters favoured by the current experimental data, in 
particular the DAMA annual modulation and XenonlOO exclusion limit. The latter 
strongly disfavours asymmetric scalar doublet DM of mass 0{ TeV) as required by 
DM-DM oscillations, while an asymmetric vector like fermion doublet DM with mass 
around 100 GeV is a good candidate for DAMA annual modulation yet satisfying the 
constraints from XenonlOO data. 
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Figure 1. A pictorial depiction of triplet scalar A partial decay giving rise to common 
origin of asymmetries in the lepton and DM sectors. 

1. Introduction 

The existence of dark matter (DM) is supported by strong gravitational evidences, i.e. 
from galaxy rotation curves, lensing and large scale structures. This implies that the 
DM particle should be electrically neutral, massive and stable on cosmological time 
scales. However its intrinsic properties are largely unknown and lead to physics beyond 
the Standard Model (SM), based on the gauge group SU{3)c x SU{2)l x f/(l)y. 
On the contrary, its relic abundance, f^DM ~ 0.23, is well measured by the WMAP 
satellite [1]. However, the mechanism that provides its relic abundance is not yet 
established. Another issue concerning SM is the origin of tiny amount of visible matter 
in the Universe which is in the form of baryons with Q^, ~ 0.04, that could be arising 
from a baryon asymmetry ns/n^ ~ 6.15 x 10^^°, as established by WMAP combined 
with the big-bang nucleosynthesis (BBN) measurements. 

In the standard cosmological picture the early Universe has gone through a period 
of inflation and then reheated to a temperature at least larger than the epoch of 
BBN. Therefore the observed baryon asymmetry of the Universe (BAU) and the DM 
component must have been generated in the thermal bath after reheating. If the 
reheating temperature is less than electroweak (EW) scale (0(100) GeV), then it is 
difficult to generate both DM and BAU [2]. On the other hand, if the reheating 
temperature is larger than the EW scale, then a handful of mechanisms are available 
which can give rise to required BAU, while leaving a large temperature window for 
creating ^dm from the DM species, which is set by freeze-out, a rather independent 
mechanism with respect to the dynamics of generating BAU. Indeed in most of the 
cosmological model, the energy density of baryons and of DM are independently 
determined. 

The fact that the energy density of DM is about a factor of 5 with respect to 
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the baryonic one could be a hint that both sectors share a common origin and the 
present rehc density of WIMP is also generated by an asymmetry. Over the years a 
large number of possibilities for asymmetric DM have been proposed [3-38]. The most 
stringent constraint on asymmetric DM candidate comes from neutron stars and white 
dwarfs in globular cluster, which exclude asymmetric scalar DM below 16 GeV [39,40]. 

In this paper we consider a relatively heavy asymmetric scalar doublet DM (SDDM) 
whose stability is provided by a remnant Z2 flavour symmetry inspired by the Inert 
Doublet Model [41-45]. Indeed the asymmetry in this model is rather natural: in 
the limit in which the number violating couphng of DM to Higgs goes to zero the Z2 
symmetry that protect the DM is elevated to a global f/(l)pQ Peccei-Quinn symmetry. 
We show that the observed relic abundance of SDDM and BAU originate naturally in 
a type II seesaw scenario [46-51], as pictorially depicted in figure 1. To accomplish the 
unification of asymmetries we extended the SM to include a SU {2)l scalar triplet A and 
an inert scalar doublet x- The partial decay width: A LL, where L is the SU{2)l 
lepton doublet, and A XX then induce the asymmetry simultaneously in both sectors. 
The lepton asymmetry is then transferred to a baryon asymmetry through the sphaleron 
transitions. In the low energy effective theory the induced vacuum expectation value 
(vev) of the same Higgs triplet gives rise to sub-eV Majorana masses, as required by 
oscillations experiments [52-54] , to the three active neutrinos through the lepton number 
violating interaction ALL + A'^HH, where H is the SM Higgs. Thus a triple unification 
of neutrino mass, asymmetric DM and baryon asymmetry of the Universe is achieved 
in a type-II seesaw scenario. 

We show that in case of a SDDM xo, the neutral component of the inert scalar 
doublet X, the asymmetry in the DM sector gets washed out below EW phase transition 
by fast oscillations between xo ^-nd its complex conjugate field xo- This sets a limit of 
the mass scale of Xo to be M^^ > 2 TeV, so that the DM freezes out before oscillations 
begin to occur. Such heavy asymmetric SDDM are quite natural to explain positron 
anomalies at PAMELA and FermiLAT, while suppressing non-observation of antiproton 
fiuxes [55]. The small number violating quartic coupling A5 of x (sec section 2.2) 
to the SM Higgs naturally provides 0{keV) mass splitting between DM particle and 
its excited state leading to inelastic interaction of detector nuclei and DM. Indeed by 
definition during inelastic scattering a DM particle that scatters off a nucleus produces 
a heavier state. It has been introduced by [56] for reconciling the annual modulation at 
DAMA [57] experiment and the null results at other experiments, e.g. [58-66]. Here we 
re-investigate [67] the compatibility of the SDDM, explaining the DAMA signal with 
the most recent exclusion bounds of CDMS-II [68], CRESST-II [69] and XenonlOO [70]. 
The SDDM appears to be strongly constrained by the new XenonlOO data. 

In analogy to SDDM, we discuss a similar model where the DM candidate is 
given by a vector like fermionic doublet, odd under a Z2 fiavour symmetry, with mass 
O(IOO) GeV, and hence to be caUed fcrmion doublet DM (FDDM). It will arise that 
the asymmetric inelastic FDDM is appropriate to explain the high precision annual 
modulation at DAMA while satisfying the latest constraint from XenonlOO experiment. 
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In that case the small mass splitting arises through a small Majorana mass of the dark 
fermion doublet given by the triplet A. 

The outline of the paper is as follows. The next section presents the particle physics 
content of the scalar triplet model which achieve a triple unification of neutrino mass, 
asymmetric inert doublet (scalar and fermion) DM and the observed BAU. After briefly 
commenting about the neutrino mass we describe the most general scalar potential 
for triplet scalar, SDDM and SM Higgs. We also address the issue of generating 
asymmetries in case of a vector like FDDM model. In section 3 we broadly discuss 
the constraints on asymmetric doublet scalar and fermion dark matter. In section 4 
the asymmetries in the baryonic and DM sector are computed and possible wash- 
out mechanisms are discussed. The Boltzmann equations are solved numerically using 
Monte-Carlo-Markov-Chain (MCMC) techniques with CP asymmetries and branching 
fractions as free parameters of the theory. Section 5 presents the constraints for inelastic 
scattering on the model parameter space from the current direct search experiments, 
using bayesian inference and marginalising over the velocity distribution of the DM 
particles. We then come to the concluding remarks in section 6. The technical details 
about bayesian analysis and MCMC are left to Appendix A. In Appendix B we show 
the results of the analysis of inelastic scattering in the case of the standard maxwellian 
halo and fixed astrophysical parameters. 



2. Scalar Triplet Model providing dcirko-lepto-genesis with Non-zero 
Neutrino Masses 

It is known that the bilinear L- violating couphng (AL = 2) of scalar triplet to lepton 
and Higgs leads to neutrino mass via type II seesaw. Moreover, the out-of-equilibrium 
decay of triplets through the same coupling can give rise to lepton asymmetry in the 
early Universe [71,72]. Here the additional decay of scalar triplets to a pair of inert 
scalar doublets (%) or a pair of vector like inert fermion doublets (ip) simultaneously 
explain the asymmetries in visible and dark sector. In our convention the scalar triplet 
is defined as A = (A++, A+, A°), with hypercharge Y = 1. 

2.1. Triplet Seesaw and Non-zero Neutrino Masses 

Since the lepton number is a conserved quantum number within the SM, the masses of 
neutrinos arc exactly zero upto all orders in perturbation theory. On the other hand, 
oscillation experiments confirm that the neutrinos are massive, however small, and hence 
they mix among themselves. This non-trivial result can be minimally explained by 
incorporating a heavy triplet scalar A to the SM of particle physics. The lepton number 
violating (AL = 2) interaction of A with SM fields is given by the Lagrangian: 

L D MlAtA + ^ [iiH^^HH -h fafsAL^L^ + h.c] , (2.1) 
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where H and L are the SM Higgs and lepton doublets respectively. After the EW phase 
transition A acquires a small induced vev, given by 

= -'"vm- 

where v — {H) — 246 GeV. The vev of A is required to satisfy 
_ M^y _ 1 + 2a;2 
^= M| cos^^^TTl^^^' ^^-^^ 
where x = {A) /v. The above constraint implies that (A) < 0{1) GeV. The trilinear 
coupling ALL then give rise to Majorana mass matrix for three flavours of light neutrinos 
as: 

{MX, = V2M^) = . (2.4) 

Hence for (A) < 0{1) GcV a wide range of allowed values of fa/s gives rise to required 
neutrino masses. For fap ~ 0{1), the required value of (A) satisfying neutrino masses 
can be obtained by choosing ixh ~ ~ 10^^ GeV. This implies that the scale of lepton 
number violation is very high. However, in presence of an extra scalar triplet the lepton 
number violating scale can be brought down to TeV scales without finetuning, so that 
its dilepton signatures can be studied at LHC [73-75] . 

2.2. Inelastic SDDM in Triplet Seesaw Model 

We now extend the Lagrangian (2.1) by including a Inert scalar doublet x = (x^X°)^ 
and impose a Z2 symmetry under which x is odd while all other fields are even. As a 
result X does not couple to SM fcrmions and hence serve as a candidate of DM. The 
interactions between A, x and H can be given by the scalar potential: 

V{A, H, x) = Mi At A + AA(At A)2 + M^jH^fH + XniH'^Hf 

+ Mlx^x + Kix'xf + [i^hA^HH + li^A^xX + h.c." 

+ MH\'\x? + MH^X? + y [{H^xf + h.c] , (2.5) 

where we have neglected the quartic terms involving A — if — x as those are not relevant 
for our discussion since the vev of A is small. The vacuum stability of the potential 



requires Aa, Ai^, A^,. > and Al = A3 + A4 — IA5I > —2^\^\h- We further assume that 
> 0, so that X does not develop any vev. This is required in order to distinguish 
the visible matter from DM given by the neutral component of the doublet x- Hence 
the true vacuum of the potential is given by: 

{H)=v; (x) = and {A) = u . (2.6) 

Since A is heavy, its vev is small as demonstrated by Eq. 2.2 and hence does not play 
any role in the low energy dynamics. Therefore, in what follows we neglect the dynamics 
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of A in low energy phenomena. The perturbative expansion of the fields around the 
minimum is: 

and X=(^) • (2.7) 

Thus the low energy spectrum of the theory constitutes two charged scalars x^, two 
real scalars h, S and a pseudo scalar A, whose masses are given by: 

r2 



=2\hv\ 

Ml =M2+(A3 + A4 + A5)" 



2 ' 

Ml =M2 + (A3 + A4-A5)y. (2.8) 

Depending on the sign of A5, either S 01 A constitutes the DM. Let us assume that A5 is 
negative and hence S is the lightest scalar particle. The next to lightest scalar particle 
is then A. The mass splitting between them is given by 

AM^ = M| - Ml = Agv^ . (2.9) 

Prom which we can deduce the coupling 

A.^^, (2.10) 

V 

where 5 = Ms — M^. This plays a key role in the direct searches of DM as we will 
discuss later. Notice that in the limit A5 — )■ in the scalar potential (2.5), there is no 
mass splitting between S and A and the two degrees of freedom can be re-expressed as 
Xo and its complex conjugate Xo- In this limit we discuss the asymmetry between xo 
and Xo via the decay of the triplet A. We then derive upper bound on DM number 
violating processes, namely XX ~^ H^H^ involving the coupling A5. The smallness of A5 
can be attributed to the breaking of a global U{1)l symmetry under which x carries a 
lepton number +1. Indeed in the absence of ^h^^HH and X^{H'^ x}^ -\-h..c. terms in the 
Lagrangian, there is a global U{1)l symmetry. The other parameters and 11^, which 
involves in the DM number violating processes xX — A — )■ HH and xX H ^ ff, f 
being the SM fermion, are not necessarily to be small as these processes are suppressed 
by the large mass scale of A. 

2.3. Inelastic FDDM in Triplet Seesaw Model 

Let us replace the inert scalar doublet x by a vector like fermion doublet i/j = {i/jdm,'^-) 
of hypercharge Y = —1/2. With the same Z2 symmetry, under which ip is odd, the 
neutral component of ip i.e. ipuM can be a candidate of DM. The relevant Lagrangian 
including the triplet scalar A is: 

-C D M^AtA + Mn^i/j + ^ [i^hA^HH + fapAL^Lf, + g^^^ + h.c] , (2.11) 
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where Md ~ C(IOO) GeV is the Dirac mass of ip. The bihnear DM couphng Aipip can 
be re-expressed as: 



-2^ 



^V2(V'i^-A++) + (V'i^DM + V'&M^-) A+ - V2(^&mV'dmA°)] , 

(2.12) 

where we have used the matrix form of the triplet scalar: 

A++\ 

_ • (2.13) 

After EW symmetry breaking the neutral component of A acquires an induced vev and 
hence give rise a small Majorana mass to ip, m — V2g{A'^). 

Therefore the Dirac spinor ■0dm can be written as sum of two Majorana spinors 
(V'dm)l and (■0DM)ii- The Lagrangian for the DM mass becomes: 

— .^DMmass = Md 




(^DM)L(^DM)i? + (^DM)i?(^DM)L 



+ m 



(V^Dm)l(^Dm)l + (^DM)ij(0DM)i?J • (2.14) 

This implies there is a 2 x 2 mass matrix for the DM in the basis {(V'dm)l) (V^dm)_r}- By 
diagonalising it two mass cigcnstatcs {iPy)m)i and (V'dm)2 arise, with masses M^^ = 
Md — m and = Md + m. Thus the mass difference between the two states 

6 — 2m ~ C(IOO) keV is required by the direct search experiments. We will come 
back to this issue while discussing inelastic scattering of DM with nucleons. Note that 
in this case the inelastic scattering of DM with nucleons ( i.e. (■0dm) i-^ — >■ (■0dm)2-^) 
via SM Z-exchange dominates to elastic scattering, as in the case of the scalar candidate. 

Now we will briefly comment about the dark matter asymmetry. Similar to the 
decay of A — > xXt the decay of A — )■ ipip will produce an asymmetry in the dark sector. 
Since ip is odd under a Z2 flavor symmetry, it will not couple to any other SM fields 
and hence the asymmetry will remain in ■^dm for ever, namely in this case there are no 
strong wash-out processes. 



2.4. Sub-eV Neutrino Mass versus keV Majorana mass of FDDM 

Notice that in case of FDDM, the induced vev of A introduces two mass scales. One 
is the Majorana mass of neutrinos, i.e. Mj, = \/2f{A^) ~ 0{1) eV and other is the 
Majorana mass of DM, i.e. m = \f2g{A^) ~ C(IOO) keV. This implies a hierarchy 
between the two couplings / (third term in Eq. 2.11) and g (fourth term in Eq. 2.11) 
of the order of O(IO^) in order to explain the triple unification of neutrino mass, 
asymmetric DM and BAU. 
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3. Constraints on Asymmetric Inert Doublet (Scalar and Fermion) DM 

3.1. Constraints on SDDM from Oscillation 

In case of SDDM the two states Xo and its complex conjugate Xo can be written in terms 
of the mass eigenstates S and A: 

\Xo) = ^{S + iA), 

\Xo) ^ ^{S - zA) . (3.1) 
The state |xo) at any space-time point {x,t) is given by 

\(j){x, t)) = ^ L-'iEst-ksx)^g^ ^ .^+iiE^t-kAx) ^ (3 2) 

where Es = \/ks + ^1 aiid Ea = \Jk\ + M\ arc the energy of S and A respectively. 
The probability of \xq) oscillating into \xo) is then given by 

^lx(.>^lx(.> = \ {Xo\^{^,t))? ■ (3.3) 
Using Eqs. 3.1 and 3.2 the probability of oscillation takes the form: 



Pi 



lxo>->-lxo> 



1 r 

4 



2 - Q-i[{Es-EA)t-{kA-ks)x\ _ ^+i[{Es-EA)t-{kA-ks)x\\ _ (3 4) 



Above the EW phase transition there is no mass splitting between the two mass 
eigenstates -S" and A, therefore Ms = Ma, Eg = Ea and ks = kA- As a result from 
Eq. 3.4 the probability of oscillation is null: 

^lxo>^lxo> = • (3-5) 

Below the EW phase transition the DM number violating term ^ {^{Wx)'^ + h.c.^ 
produce a mass splitting between the two mass eigenstates S and A. From Eq. 3.4 the 
probability of oscillation becomes: 

1 



Pi 



lxo)->|xo> 



2 



/AM2(t-tEw)\ 

COS 7— 



(3.6) 



where we have assumed Es ~ Ea ~ E, which is a good approximation for a small 
mass sphtting. In the following we will consider a mass sphtting of C(keV), which 
imphes A5 ~ 10~^. We also normalise the time of evolution from the time of EW phase 
transition, so that t — t-Ew, P\xo)-^\xo) — 0- Below EW phase transition the time of 
oscillation from xo to Xo can be estimated as 

In the relativistic limit the energy of the DM particle ~ T, where T is the temperature 
of the thermal bath. Hence the oscillation time can be given as: 
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X 



A 



Figure 2. Tree level and self energy correction diagrams for the production of CP 

asymmetry in leptogenesis. 



On the other hand, in the non-relativistic hmit the energy of the DM particle E ~ M5. 
Thus for Ms ~ 100 GeV, the time of oscillation is again similar to relativistic case. This 
implies that xo oscillates rapidly to xo- In this case if xo is in thermal equilibrium then 
during each oscillation there is a leakage of asymmetry through the annihilation channel 
XoXo — >■ SM particles. Alternatively to keep the generated asymmetry intact xo should 
freeze-out before it oscillate to Xo- In other words, the mass of Xo should be given by 



where ~ 20, which determines the epoch of freeze-out. From the above equation we 
see that to get an asymmetric SDDM one should have M^^ > 2 TeV. 

3.2. Constraints from Collider 

Since the DM (scalar or fermion) is a doublet under the SM gauge group, it couples to 
the Z boson. As a result they can change the invisible decay width of the later unless 
the mass of DM is greater than half of Z-boson mass. This gives a lower bound on the 
mass scale of either SDDM or FDDM to be > 45 GeV [76]. 

4. Developing asymmetries in the lepton and DM sectors 

If the triplet A is heavy enough as required by the seesaw, then it can go out-of- 
equilibrium even if the gauge couplings are [71,72,77,78]. In such a case the 

out-of-equilibrium decays of A — > LL and A ^ XX (V'V') produce asymmetries in 
visible and dark sectors respectively. The CP asymmetry for the two sectors arise 
via the interference of tree level decay and self-energy correction diagrams as shown 
in figure 2 and 3 respectively, for the scalar DM case, but totally analogous for the 
fermionic doublet. 

Considering the inert scalar doublet as a reference for the scalar potential, from 
Figs. 2 and 3 we see that the CP asymmetry requires at least two triplet scalars. Hence 
in presence of their interactions, the diagonal mass in Eq. 2.5 is replaced by: 



M^^ ^ x/Tew, 



(3.9) 




(4.1) 
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H 



Figure 3. Tree level and self energy correction diagrams for the production of CP 
asymmetry in generating asymmetric DM. 



and the trilinear couplings hh^^HH + n^A'^xX + h.c. in the scalar potential (2.5) 
become: 



a=l,2 



In Eq. 4.1, the mass matrix is given by: 



M 



± 




—iCi2 
M| - iG 



22 



where 




(4.2) 



(4.3) 



(4.4) 



with C, 



ab 



r:,Mb, and Ca 



^Ma- Solving the mass matrix 4.3 one gets two mass 
eigenstates C,i_2 = + with masses Mi and M2. The complex conjugate 

of ^^''2 given by .^f2 = ^12^1 + ^12^2- Note that C,^ and states are not CP 
eigenstates and hence their decay can give rise to CP asymmetry Wc assume that there 
is no asymmetry, either in the visible sector or in the dark sector, at a temperature 
above the mass scale of the triplets. The asymmetries are generated in a thermal bath 
by the decay of these triplets. If we further assume that the mass of is much less 
than the mass of ^2 then the final asymmetries in visible and dark sectors will be given 
by the decay of as: 

ei = 2 [Bt{^- U) - Br(^+ ^ 

= 2 [Br(er ^ XoXo) - Bri^t ^ X*oXo) 
where the front factor 2 takes into account of two similar particles are produced per 
decay. From Figs. 2 and 3, the asymmetries are estimated to be: 



= ev 



— ^dark ) 



(4.5) 



Im ( liixlJ'*2x 



J2al3 flaf3f2a/3 



and 



Im (, 



1^1x1^2), 



87r2 (M| -M2) 

flafif2a/S 



Ml 



87r2 (M| - Ml 



Ml 



Pi 



(4.6) 



(4.7) 
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where Fi = Fn. 

In a thermal bath these asymmetries evolve as the Universe expands and settle to 
a final value as soon the relevant processes go out of equilibrium, i.e. 

Fi = ni{ai\v\) < H{T) , (4.8) 

where H{T) is the Hubble scale of expansion. As a result the Yields in both sectors can 
be written as: 

s 

Ti 

^ ^ = e^X^r^^ , (4.9) 

where — n^-/s = n^+/s, s — 2{7r^ /A5)g*T^ is the entropy density and r)L,r)^ are 
the efficiency factors, which take into account the depletion of asymmetries due to the 
number violating processes involving x, L and H (this holds also for the fermionic inert 
doublet, hence we can replace x ~^ DM label). At a temperature above the EW phase 
transition a part of the lepton asymmetry gets converted to the baryon asymmetry via 
the SU{2)l sphaleron processes. As a result the baryon asymmetry [79] is: 

= - frtlr ^^ = -^^^^^ ' (4-10) 

14n + 9m 

where n is the number of generation and m is the number of scalar doublets, leading to 
Sum — 0.53, 0.55 for scalar DM and fermionic DM respectively. 

As introduced in Sec 2.2, in the case of the scalar doublet DM, the asymmetry may 
strongly washed out, if kinematically allowed, by the DM number violating processes: 
XX ^ A — )■ HH., XX ~^ H'^H'^ (contact annihilation through A5 coupling) and 
XX ~^ H ^ f f . The reduced cross-section for the former process is given by: 

1 l/^xlVg | 

^TT (S - Ml 

where s is the centre of mass energy for the process: XX — ?■ A ^ HH . Below the mass 
scale of the triplet this process is strongly suppressed. On the other hand in case of the 
contact annihilation of x's the reduced cross-section is given by 

As a result the reaction rate is given by F^ = {''ixl^^x^-, where the reaction density is 

T 



a(xx^A^M) = ^l'^>^^, (4.11) 



7x 




ds^sKA^]ax, (4.13) 



647r4 

and the equilibrium number density of % is 

„y.-|^,,(^), (,14) 

where g^oi is the internal degrees of freedom and s is the usual Mandelstam variable 
for the center of mass energy. In Eqs. 4.13 and 4.14, Ki and K2 are modified Bessel 
functions. In fig. 4 we compare the rate of the process xX ~^ H^H^ with the Hubble 
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Figure 4. The scattering rate of the process XX ~^ for different values of A5 

is compared with the Hubble expansion rate. For illustration purpose we have used 
= 2 TeV and Mi = IQi" GeV. 



rate by taking three values of A5. We see that for A5 < 10^^ (blue dotted line), the 
scattering rate remains out-of-equilibrium through out the epoch. For larger values of A5 
the scattering process comes to equilibrium at late epoch. At around z = Mi/T ^ 10'', 
which implies M^/T = {M^/Mi)z ~ 1, the scattering rate of the process sharply drops 
as it is expected and does not depend on the value of A5. This argument also holds 
in case of the process xX H ^ f f ■ However, the rate of the scattering is further 
suppressed by the mass scale of Higgs. 

Returning to the general case of doublet DM, both scalar and fermion, from Eqs. 4.9 
and 4.10 the DM to baryon ratio is given by: 

7^ = ^ ' ^^-^^i 

\Ib odm mp eL Vl 

where ~ 1 GeV is the proton mass. From this equation it is clear that the dependence 

of Vbm/vl on the mass of DM goes as l/moM- Hence for a C(IOO) GeV scale DM, the 

required efficiency factor for DM is two orders of magnitude less than the case of lepton 

provided that the CP asymmetries are equal on both sectors (see the end of section 4.2). 

For example, from Eqs. 4.6 and 4.7 we notice that the CP asymmetries are identically 

equal if: 

^ = E/iW.*«,. (4.16) 

1 a,/S 
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In what follows, section 4.2, we solve numerically the relevant Boltzmann equations 
for quasi equilibrium evolution of triplet scalars, presented in section 4.1, to show that 
the parameter space of the theory fulfills the criteria Qbm ~ 5 Ob and the observed 
BAU, while allowing a broad range of asymmetric DM masses. 

4.1. Boltzmann equations for quasi- equilibrium evolution of triplet scalars 

If the triplet {^f) decay occurs in a quasi-equilibrium state then the detailed of rj^ and 
r]L depends on the dynamics of the processes occuring in the thermal bath and can be 
obtained by solving the relevant Boltzmann equations [80, 81]. In our case the additional 
decay channel of the scalar triplet into DM particles is included. 

At first the number density of particles changes due to their decay {^f — > 
LL, HH, XX or V'V') and gauge annihilation {^^^t ^ LL, WH, x^'xi^^): W^W^, B^B^), 
where and B^ are the SU{2)l and U{1)y gauge bosons respectively. If we assume 
that the masses of the components of the triplet are same before EW symmetry 
breaking then it is fairly general to use the dimensionless variables z = Mi/T and 
= n^-/s = n^+/s. The Boltzmann equation for the evolution of density is then 
given by: 

~dr - ~Jh{z) V^«^ -^^^)~ Jh{z) \ Xl^ ) ' (^-^^^ 

where 

r„ . r.|lW , r, ^ ^ and HU) ^ ^^I^ , (4.18) 

and 7a's arc the scattering densities, described in Eqs. 4.20. The temperature 
independent decay rate of .^1 can be written as a function of the neutrino mass: 

r, = -1 1"^-'^^' (4 19) 

8n{H)^,/BZB^' ^^-'^^ 
where B^ and Bh are the branching fractions in the decay channels: ^1 — > LL and ^1 — > 
HH. Note that we have re-expressed the total decay rate ^i{fDM: fn, fz. Mi), where 
/dm = l^x/Mi for SDDM and g for FDDM and fn = I^h/Mi, as ri{m^, Bl, Bh, Mi). 
In the following we set my = 0.05 eV and therefore the total decay rate depends only 
on three variables, namely Bl, Bh and Mi. This makes a crucial decision in setting up 
the final asymmetry as we will show in section 4.2. 

For the gauge annihilation processes, the scattering densities are given by: 

licuT ^ //) = "^'^''ft'"^^ r dxv^Mzv^y , 

^ H'H) = £ dxV^Ki{zV^)r' , 

^ X'X) = ^'^^l^/f^^ /I dxV^Ki{zV^)r' , 
7(ei+er ^ W^W) = -if / dxV^Kiiz^) 



24 / 1 -I- r ' 

r(5 + 34/a;) - — (x- l)ln ' 
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7fe+^r ^ BB) 



1287r5^ 

X 



dx\/xKi{z\/x) 



r(l + 4/x) - 4(^-2) In 



(4.20) 

jj~ \ J- — / / J 

where r = y'l — 4/a; and a; = s/M^. In case of FDDM, the process corresponding to 
iiii X^X is given by: 

Mt{Qgt + bg^) foo 



1287r5^ 



/ (ia;V^is:i(zV^)r^ (4.21) 



[n 



Since are charged particles there is an evolution of the asymmetry: 1^^ 



?r ~ "^(.t^l^ decay and inverse decay of particles. The evolution of 

(4.22) 



1^^ is described by the Boltzmann equation: 



dz 



zH{z) y zH{z) ^ ^' 



where Yj = {rij — nj)/s, with j = L,H,x {"ip) and 



-r «i 

>- ID — D -^eq 



and = , 

J- 1 



(4.23) 



where Xj 
equation: 



Uj/s. The evolution of the asymmetries Yj is given by the Boltzmann 



dY^ 
dz 



2{ 



D 



zH(z) 



D 



zH{z) 



_3^2y' 

zH{zr\ 



x: 



eq 



^ zhU xr^^^y 



(4.24) 



where Ts = •js/nl- is the scattering rate involving the number violating processes, such 
as XX ^ C ^ HH, LL — )■ { ^ HH. The front factor in Eq. 4.24 takes into account of 
the two similar particles produced in each decay. 

Solving the Boltzmann Eqs. 4.17, 4.22 and 4.24 we can get the lepton (Yj,) and 
dark matter (1dm) asymmetries. Note that because of the conservation of hypercharge 
the Boltzman equations 4.17, 4.22 and 4.24 satisfy the relation: 2Y^ + Yj = 0. This 
implies: 



Y, 



3 



(4.25) 



We will follow a phenomenological approach and calculate the ratio of efficiency 
factors r/DM and r]^ (and hence also the individual efficiency) solving the set of coupled 
equations 4.17, 4.22 and 4.24. As usual the efficiency factor for the species i = L,H, DM 
is defined as: 

Vi = ■ (4.26) 

Xc 

^ T»Mi 

The free parameters of the model are the CP asymmetries for all the species, 
the dark matter mass and the triplet mass Mi. However in the remaining of the 
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Figure 5. Le/f; Likelihood function £ratio given by the ratio distribution, as a function 
of r = f^DM/^b- Right: Logarithm of the ratio hkelihood (red hne) and logarithm of a 
gaussian distribution (dashed gray curve). Both curve have the same variance. 

paper we will focus on heavy triplet, Mi ~ 10^° GeV, which has been shown to lead to 
successful leptogenesis [80, 81] for a wide range of CP asymmetries and branching ratios. 
In addition the following constraints apply: 

II 9 = ' H = 1 and le^l <2 Bj. (4.27) 

j j 

The first and third conditions ensure that all amplitudes are physical and the total 
amount of CP violation can not exceed 100% in each channel, while the second condition 
simply demands unitarity of the model. The number of free parameters therefore drops 
to 5, which we choose to be: e^, eDM, B^, -Bdm and rriDM- The numerical procedure and 
the results are described in the next section. 

4-2. Numerical solutions of the Boltzmann equations 

In principle for unlimited computational power one could use a griding method to explore 
the whole parameter space and to localise the hypervolume that satisfies Eq. 4.15. For 
the problem under scrutiny we have a 5-dimensional space: MCMC technique is well 
suited in this case since it numerically scales linearly with the number of dimension 
instead of exponentially. We follow the approach presented in [82] |. 

Defining a probability measure over the full parameter space allows to use Bayesian 
inference to assess the posterior probability distribution of all the parameters to get 
asymmetry in the dark sector as well as in the lepton sector, which gives rise to observed 
BAU. The details about bayesian statistical methods and the implementation of the 

I MCMC technique has been applied for scanning hybrid inflationary parameter space [82] and for 
investigating the multi-dimensional parameter space of supersymmetric theories, e.g. [83]. 
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Table 1. MCMC parameters and priors for the CP asymmetries, branching ratios and 
moM- All priors are uniform over the indicated range. 



MCMC parameter 


Prior 


log(mDM/GeV) 


0^5 


log(eL) 


-9^0 


log(eDM) 


-9^0 


log(i?L) 


-5^0 


log(5DM) 


-5^0 



MCMC are given in Appendix A, while for detailed reviews we refer to e.g. [84, 85]. The 
important point that has to be underhned is that through Bayes' theorem: 

P(^|X)d^ oc £(X|^) ■7r(^)d^, (4.28) 

the posterior probability density function (pdf) V{6\X) for the model parameter 6, given 
the data X, is proportional to the hkehhood of the experiment times the prior behef in 
our model 7r{9) and is sampled directly by the MCMC elements. 

The hkelihood function jC{X\9) denotes the probability of the data X given 
some theoretical prediction 9 and plays a central role in Bayesian inference. Both 
the abundances of dark matter Qdm and baryonic matter il^ a-re variables normally 
distributed around their mean values given by WMAP measurements [1]: JIdm = 
r2i:)M±o"DM = 0.227± 0.014 and = Qh±aiy = 0.0456 ±0.0016. Since we are interested 
in the ratio of the dark to baryonic matter, Eq. 4.15, we define the likelihood as the 
probability distribution of the model parameter to satisfy that ratio. The likelihood 
is therefore well described by the so-called ratio distribution [86], which is constructed 
as the distribution of the ratio of variables normally distributed with non zero mean. 
Calling r = ^Du/^b = finT-DU, Bi), the likelihood reads: 

' &(^)c(r) ^/6(r)x ^ exp ( - l/2(ng,M/^DM + ^feK)) 



r 

with: 



a(r) 



1 



2 ' 



"^DM '^b 



= -2—^ + -2 ' 



c r 



^{u) = Erf (^) . (4.30) 

with Erf being the error function. The shape of the likelihood function is depicted in 
figure 5, where it can be seen that the peak is at around r ~ 5 as it is expected. In 
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Posterior pdl 



Figure 6. Left: ID posterior pdf for the DM mass ttidm- Right: 2D credible regions 
at 68% and 95% C.L. in the {todm, '7DM/?7L}-plane. The vertical dot-dashed blue line 
denotes the bound from LEPII, while the red dashed vertical line marks the bound 
from xo ~ Xo oscillations for SDDM. All of other parameters in each plane have been 
marginalized over. 



the right plot we show that the ratio distribution is shghtly skewed with respect to a 
gaussian distribution (gray dashed hne) with the same variance. 

In addition the baryon asymmetry should satisfy the constraints from WMAP: 



r^feia^b = (6.15 ±0.25) X 10"'% (4.31 



where rib = 7.02 x SbmYl, and the density of photon and baryons are computed at 
present time. The baryon asymmetry is described by a gaussian distribution: 



£l oc exp ( — . (4.32) 



Summarizing, the logarithm of the total likelihood is given by the sum of Eqs. 4.29 
and 4.32: 

In -^asym In .^ratio 

+ ln£L. (4.33) 

MCMC techniques require a prior assumption on the probability distribution of 
the parameters, namely on vr(6'). In the absence of theoretical constraints on the model 
parameter there are a priori no constraints on these quantities. We therefore focus on 
the regions singled out by successful triplet leptogenesis in [80, 81], that is values of the 
CP asymmetries ranging from 10"^ — 1 and branching ratios from 1 to 10"^. In order 
not to support a particular scale we choose flat prior on the log distribution for each 
parameters, as described in table 1. The dark matter mass is let free to vary between 1 
GeV up to 10 TeV (even though masses below 50 GeV are excluded by LEPII). As it is 
known, if the data are not informative enough a dependence on the prior choice is left 
in the posterior pdf. 

We present in the following the results of the bayesian inference for the case of the 
asymmetric inert scalar doublet dark matter (SDDM). 
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Figure 7. Left: 2D posterior pdf in the {BuM,VDM/'r]L}-p^3.ne. Right: Same as left 
in tlie {Sl, 'yDM/'yij-plane. The credible regions are given at 68% and 95% C.L.. All 
of other parameters in each plane have been marginalized over. 



In the left panel of figure 6 we show the ID posterior pdf for moM, while all other 
parameters are marginalized over. We see that all the mass range from 1 GeV up to ~ 4 
TeV can lead to successful leptogenesis, namely Yl ~ 10^^" and an asymmetric dark 
matter satisfying the ratio Eq. 4.15. The vertical blue line denotes the bound from the 
Z decay width and masses below 50 GeV are excluded, while the region on the left of 
the red vertical line is excluded by Xo~Xo oscillations. As one can see from the posterior 
pdf the most favoured region is at mDM ~ 10 GeV, while there are candidates at 100 
GeV with smaller statistical significance. With even less probability but still viable 
are candidates at TeV, which is the range of interest for SDDM. On the right panel 
of figure 6 the 68% and 95% credible region are shown in the {mDM, "/^DM/^^Lj-plane. 
From there we see that for DM mass up to around 500 GeV, the preferred values of 
the ratio //dm/^/l remains constant to be around 10-50 as these are easily compensated 
by the small CP asymmetry ratio eBu/^L- However, for DM masses above 0{TeV), 
cdm/cl is not sufficiently small to compensate with large ?7dm/'7l and therefore, the 
preferred values of tj-dm/vl remains to be around unity. Alternatively for a given DM 
mass, smaller values of tidm/til are allowed at 95% C.L. for e-Du/^L > 1- 

The above preferred values of the parameter space can be understood from the 
Boltzmann equation 4.24 in which the input parameters are the CP asymmetries and the 
branching ratios. Having chosen a mass of the scalar triplet of 10^*^ GeV the dominant 
processes that regulate the boltzmann equations are the decay and inverse decay, and 
fundamental quantities are the branching ratios. In figure 7 we show the correlation 
of rjDu/VL versus -Bdm and Bl respectively in the left and right panels, within the 
68% and 95 % credible regions. We see that large efficiency ratio ?7dm/'7l is preferred 
when Bl ^ 1 and small -Bdm — ^ 10~^. This is because larger the value of Bl (which 
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Figure 8. Left: 2D posterior pdf in tlie {j/dm, ?/L}-plane. Central: 2D posterior pdf 
in tlie {i^DM, 'yDMj-plane. Right: same as left in the {Bl, ?7L}-plane. The credible 
regions are given at 68% and 90% C.L.. All of other parameters in each plane have 
been marginalized over. 



implies smaller is the -Bdm as = 1 with i=L, H ,DM) the larger is the washout 

due to inverse decay and hence leads to small rjL. On the contrary smaller is the -Bdm 
the washout effect is small due to inverse decay and hence large //dm- Note that in 
either case the production of asymmetry is proportional to Fi oc 1/ ^BiBh. Therefore 
when Bl approaches towards 10~^ the asymmetry {Yl) as well as the efficiency {til) get 
increased. On the other hand when -B^m approaches towards 1, which implies small -B^, 
the asymmetry Fdm gets increased but efficiency gets decreased. These behaviours of 
?7dm and rji can be confirmed from figure 8 where we have shown the 2D credible regions 
at 68% and 95% C.L.. The extreme left one, which constitutes the summary of middle 
and right ones, reveals that a successful asymmetric dark matter and lepton asymmetry 
can be generated with small rji and large r/oM- In other words, large Bi and small -Bdm 
are required in favour of the observed BAU and asymmetric dark matter. 

For sake of reference we report the preferred values of the input CP asymmetries. 
The preferred values range between 10~^— 10~^ for eoM, respectively to -Bdm = 10~^— 0.5. 
A more tighter range is selected in the case of the lepton CP asymmetry: 10~^ — 10~^ 
again for Bi ranging from its extremal values. We remark in addition that for large 
-Bdm and small Bl and masses around 50 GeV, the CP asymmetry in the DM sector 
can be be larger by an order of magnitude with respect to to compensate the small 
value of rjBu/VL- 

Regarding the inert fermionic doublet DM (FDDM) candidate, the discussion is 
very similar. We have verified that there are no substantial differences in the selected 
ID and 2D credible regions as the Boltzmann equations are same in both cases. A small 
difference comes from the internal degrees of freedom which makes the equilibrium value 
of a FDDM different from a SDDM. Therefore, in case of FDDM, the allowed mass range 
goes up to a few TeV starting from 50 GeV as shown in figure 6. 

In figure 9 we show the behavior of the Yields for leptons, Higgs, DM, scalar triplet 
and for two particular points. The first point in the parameter space is shown in 
the left panel, which leads to a successful model for FDDM with a mass of ~ 86 GeV, 
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Figure 9. Left: Absolute value for the Yield of leptons (cyan solid), DM (dotted 
magenta), Higgs (dashed black), ^ asymmetry (solid red) plus scalar triplet abundancy 
(black solid), for a successful point with ttidm — 86 GeV, -B^ = 0.09, i?DM = 4.1 x lO""*, 
cl = 2.6xl0-^ EDM = 1.1x10-8 which leads to r = f^DMM = 4.75, Yl = 1.7x10-1° 
and ?7dm/7?l — 6.53. Right: Same as left for todm — 2 TeV, Bj^ = 9.5 x 10^'^, 
Bum = 2.6 x lO'^, = 7 x lO-\ eoM = 1-2 x 10^^, I^dmM = 5.4 and 
Yl = 1.6x 10^1° and tiom/iIl — 0.86. The \Yi\ are rescaled in terms of CP asymmetries. 

r ~ 4.8 and Yl = 1.7 x 10~^°. The second point in the parameter space is depicted in 
the right panel and accounts for a SDDM with m-DM ~ 2 TeV, r ~ 5.3 and successful 
baryon asymmetry, Yl = 1.6 x IQ-^*^. The details about the parameters are given in 
the caption. These two points are representative of the behavior discussed above. In 
particular, for the left panel the branching ratios are Bl = 0.09 and 5dm = 4.1 x 10~^, 
which implies small tjl and large ?7dm- Therefore, the ratio of rjDyi/rjL is maximum 
and can be confirmed from figure 7. For the figure in right panel the branching ratios 
are both small Bl = 9 x 10~^ and 5dm = 3 x 10~^, which implies t^dm and rjL are 
comparable. As a result the ratio tidm/vl ~ 0.9 and the large DM mass is compensated 
by the very small CP asymmetry ratio. This behavior can be confirmed from figure 7. 

For checking the consistencies, we have investigated the behavior of efficiency factors 
in case of equal CP asymmetries in DM and lepton channels. There are two interesting 
results that come out from the MCMC run, shown in figure 10. In the left panel the ID 
posterior pdf for the DM mass is depicted. We note that equal asymmetries lead to a 
upper bound to the asymmetric dark matter mass of C(50) GeV and can be applied to 
the case of FDDM. This is some how expected because the ratio of CP asymmetries can 
not compensate the increasing DM mass and therefore around 100 GeV there are no more 
candidate which can fulfill the DM to baryon requirement. This can be seen from the 
central panel where the 2D credible regions in the {moM, V'dm/vl}-V^^^g are shown: the 
ratio of efficiency drops very rapidly as soon the DM mass increases. As a consequence 
we found the known results of a light DM mass, around 10 GeV, being favoured. In 
the right panel the 2D credible regions are shown in the Bl — -BDM-plane. We see 
that the preferred branching ratios are similar in magnitude as expected and therefore 
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Figure 10. Left: ID posterior pdf (black solid line) for the DM mass moM- The 
vertical dot-dashed blue line denotes the bound from LEPII. Middle: 2D credible 
regions at 68% and 95% C.L. in the {mDM, ??DM/?7.L}-plane. Right: Same as middle in 
the {Bdm, Sij-plane. All of other parameters in each plane have been marginalized 
over. 

the efficiencies are comparable. In other words the preferred ratio tj^m/vl ^ 0.1 for 
rriBM < 50 GeV. 



5. Inelastic scattering and Asymmetric DM 

In this section we briefly recall the definition of event rate in a detector and discuss the 
features of inelastic scattering. For each experiments we describe the likelihood functions 
used in the data analysis and the choice of priors. In the case of direct detection signals, 
the only free parameters of the model are m-ou and the mass splitting 6. For details on 
the experimental set up and bayesian inference we refers to Appendix A and [87]. 



5.1. Experiment description, likelihoods and priors 

The direct detection experiments aim to detect or set limits on nuclear recoils arising 
from the scattering of DM particles off target nuclei. The energy En transferred during 
the collision between the incident particle with mass moM and the nucleus with mass 
is of the order of the keV for a mean DM velocity of f/c ~ 10~^ in the Galactic 
halo. The differential spectrum for such recoils, measured in events per day/kg/keV, is 
given by 

where da/dE^ encodes all the particle and nuclear physics factors, Pq = PdmI-Rq) is 
the DM density at the Sun position and ri{E{i,t) is the mean inverse velocity of the 
incoming particles that can deposit a given recoil energy Ej^: 

,(E,,t)=f dV?(fWl, (5.2) 
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where the velocity v is taken with respect to the Earth frame. The quantity I'min is the 
minimum velocity needed to lead to a recoil inside the detector: 

where /i„ the WIMP-nucleus reduced mass and Mj^ is the miclcon mass. 5 denotes the 
mass splitting between the DM particle and the excited state, therefore proportional 
to A5 in case of SDDM and Majorana mass in case of FDDM. The value required by 
preserving the asymmetry, A5 ~ 10~^ leads precisely to S ~ 100 kcV, the right order 
of magnitude for inelastic scattering in case of SDDM. On the other hand, in case of 
FDDM, 5 ~ m, the Majorana mass of FDDM. 

The total event rate per unit detector mass is obtained integrating in a given energy 
bin [Ei,E2] Eq. (5.1), 

dEn e{En) ^ 
lEi dEu 

where e is the energy dependent efficiency of the detector. The expected number of 
event observed in a detector is given by: 

S = M^,,TR{t) , (5.5) 

where M^^t is the detector mass and T is the exposure time. 

For some detectors, like scintillators, the recoiling nucleus may loose energy by 
collisions with other nuclei, hence in form of heat, or through collisions with electrons, 
which create scintillation light. The observed energy released in scintillation light 
(typically expressed in keVee) is related to the nuclei recoil energy through the quenching 
factor g, ^scint = qEr. 

The particle physics cross-section for coherent inelastic scattering S M ^ A J\f, 
mediated by the Z exchange on t-channel, is parameterised as: 
da Mm Gl 



R{t) = j^' dEn e{En) ^ , (5.4) 



dEn 2^1 2nf, 
Mj^ G 



{A-Z)U + Zfp) F\Er) 



27r + ^ ~ ^)fF\ER) , (5.6) 

where Z and A are respectively the number of protons and the atomic number of the 
element, Gp is the Fermi constant and sin^ is the Weinberg angle. For SU{2)l 
doublets with hypercharge 1/2 the couplings to proton and neutron are different, i.e. 
fn 7^ fp in contrast to the standard elastic scattering. We note that the cross-section 
is no longer a free parameter. The nuclei form factor F'^{Eb) characterizes the loss of 
coherence for non zero momentum transfer and is described as the Helm factor [88, 89]. 

Regarding ri{Eji, t) in equation 5.2, we consider the velocity distributions generated 
by a cored isothermal and the NEW [90] density profile marginalized over the 
astrophysical variables {vq, v^sc and p©), while the standard maxwellian halo is presented 
in Appendix B for sake of reference. Given a DM density profile and assuming 
equilibrium between gravitational attractive force and pressure, a corresponding velocity 
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distribution arises from the Eddigton formula [91]. We consider these two DM density 
profiles because they lead to velocity distributions that differ mainly in the selected 
mean and escape velocity values [87], which are fundamental quantities in inelastic 
scattering. Indeed the splitting factor 6 in Eq. 5.3 means that only the very high 
velocity particles will have enough energy to produce a recoil in the detector. This can 
be seen re-expressing it in terms of target and DM masses 

/v\2 2S{M^f + mpM) , 

In addition experiments with heavy nuclei will have a large sensitivity to the high tail 
of the velocity distribution. We therefore discuss the DAMA experiment because of the 
Iodine, CRESST-II on W and XenonlOO. The Germanium is more sensitive to particle 

of mass of the order of 50-70 GcV but we consider the dedicated analysis for inelastic 
scattering [68]. We do not consider the CoGeNT [92] experiment since only very light 
DM can account for its excess. We do not consider Zeplin-III [93] since it has analogous 
sensitivity than CRESST-II. 

CRESST-II CRESST is a cryogenic experiment running at the Laboratori Nazionali 
del Gran Sasso. The 33 detector modules are made by GaW04 crystal, each of a mass 
of 333 g. In this analysis we consider the second run of GRESST, carried out in 2007 
and in particular the data on Tungsten [69] . These data are obtained with two detector 
modules, leading to a total exposure of 30.6 kg days on W after cuts, in the energy range 
of 10-40 kcV. Three events have been seen, which are compatible with the expected 
background, mainly from neutrons, of ~ 0.063 kg days. We therefore use a background 
value B = 3 with as = 10% of B. The likelihood is described by the poisson probability 
of seeing three events for a given theoretical prediction S and a given background B: 

lnCcressMS,B) ^ -{S + B) + 3ln{S + B) . (5.8) 

The effective likelihood we used in the analysis is marginalized numerically over the 
background: 

ln£g,,,, = / dB ln£cresst(3|5, B) p{B) , (5.9) 

where p{B) is the probability function of the background, modeled as a gaussian 
distribution. The invariant 90^% confidence level, based on the yS-signal, corresponds 
to Ax^ ~ 3.34. 

CDMS on Germanium The GDMS collaboration has published a dedicated analysis 
for inelastic DM [68], which we use for constructing the likelihood of the experiment. 
The total exposure is 969 kg days and the energy range is from 10 keV up to 150 keV. 
From AEi = 10 — 25 keV 8 events has been found, with an expected background of 
5.88lx;75, while in the remaining energy range, AEi = 25 — 150, 3 events survive all 
the cuts and have an expected background of 0.93lo;3g. The background accounts for 
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all surface events and cosmogenic particles. For the detector efficiency we used the red 
dotted curve presented in figure 5 of [68]. 

The total likelihood is the sum of the contribution from the two energy range AEi 
and AE2. Each partial likelihood follows the poisson distribution: 

\nCcDMs{ll\S,B) = IuC^eA^IS^B) + \nCAEA^\S:B) , (5.10) 

with 

ln£ABi(8|-5,S) = -{S + B) + S\n{S + B) , 

\njC^EMS,B)= -{S + B) + 3\ii{S + B) . (5.11) 
We then marginalise numerically over the background: 

^CDMS = / CcBMs{n\S, B) p{B), (5.12) 
Jo 

to get the effective likelihood we use in computing the exclusion bound. 

The 905% confidence interval corresponds to Ax^ = 2.5, obtained considering that 
in the whole energy range there are 11 events with an expected background of 6. 

DAMA The DAMA likelihood for the modulated rate is described in [87] and follows 
a Gaussian distribution: 

ln>CDAMA = - E o 2 ' (5-13) 

where and sf^^ are the theoretical and the mean observed modulation respectively in 
the iih. energy bin, o"j is the associated uncertainty in the observed signal. We use in 
this analysis the 12-bin data from figure 9 of [94]. The quenching factors ^Na and qi are 
taken to be free parameters in our analysis, which we vary over their respective allowed 
range [95,96], as reported in table 2. In addition we require that the unmodulated 
predicted signal does not overcome the total unmodulated rate in figure 1 of [94] , namely 
in each energy bin the predicted total rate should be at most equal to the measured 
rate. 

XenonlOO As far as it concerns XenonlOO (XelOO hereafter) experiment, we use the last 
data release for inelastic scattering [70] . The likelihood is given by a Poisson distribution 
for three seen events times a gaussian which takes into account the uncertainties on the 
scintillation efficiency Leff . These latter however affect only the low mass DM region. In 
addition the uncertainties over the background are marginalized over. For details about 
this experiment we refer to [87] . 

Choice of priors and pills of hayesian inference As for the bayesian inference we follow 
closely the approach of [87]. We consider a full bayesian analysis without discussing 
profile likelihoods: indeed in case of informative data (as for DAMA experiment) the 
posterior pdf and the profile likelihood are equivalent, while in the case of exclusion 
bounds, in order to be insensitive on the choice of priors, we use the invariant bound 
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Table 2. MCMC parameters and priors for the model parameter space and 
experimental systematics (nuisance parameters). All priors are uniform over the 
indicated range. 



Experiment 


MCMC parameter 


Prior 


All 


log(mDM/GeV) 


0^5 


All 


5/keV 


^ 200 


DAMA 




0.2 ^ 0.4 


DAMA 


Qi 


0.06 0.1 


XenonlOO 




-0.01 ^ 0.18 



for the xs% credible region, as described in Appendix A, which can have a bayesian 
interpretation in terms of probability for the S signal. 

Having specified the likelihood functions for each experiments, the only missing 
element for Bayes theorem, Eq. 4.28, is the choice of priors. As in the previous section, 
the prior on the DM mass is chosen flat on a logarithmic scale on the same range (note 
that it is the only parameter in common with the asymmetry generation MCMC), while 
for 5 we chose a fiat prior, since the scale of this parameter is known by the requirement 
of inelastic scattering. The range is given in table 2, together with the priors for the 
systematic parameters in each experiment. 

Astrophysics In addition to the candidate mass moM, the mass splitting 5 between the 
DM and its excited state and the nuisance parameters in the experimental set-ups, two 
further free parameters are used to characterise the DM velocity distribution: the virial 
mass of the DM halo, and its concentration. These additional parameters are, however, 
also constrained by astrophysical observations on the velocity of the local standard at 
rest, Wo, on the escape velocity for the DM halo, v^sc and on the DM density at the 
sun position p©, which all enter in equations 5.1 and 5.2. The gaussian priors and the 
astrophysical likehhood are given in details in [87]. Only in the case of Maxwellian 
velocity distribution we do not vary the astrophysical observable in their allowed 
range but keep them fixed at their mean values, which are Vq — 230 kms~^ [97,98], 
^;esc = 544 kms-^ [99, 100] and = 0.4 GeVcm-^ [101, 102]. 

5.2. Results for scalar and fermionic candidates 

In this section we present our inference analysis for the considered experiments. Before 
coming to the results for the scalar and fermionic candidate we would like to make few 
general comments about inelastic scattering mediated by the Z boson and the DAMA 
region. 

Figure 11 shows the preferred DAMA region for inelastic scattering in the plane 
{6, m-Du} for a NFW density profile. In the left plot only the modulated signal is 
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Figure 11. Left: 2D marginal posterior for DAMA in the parameters space {<5, todm} 
and the NFW density profile. Right: Same as left with the additional constraint of the 
total unmodulated rate. 



considered, while on the right-hand plot we add the additional constraint on the total 
rate: it follows that the region with small 5 ~ 20 keV and large masses, (9(10) TeV, 
is excluded. In addition the mass range around 50 GeV and small mass splitting is 
disfavoured, while an island at lower masses and splittings survives. The same behavior 
is retrieved for the cored isothermal halo. In figure 12 the dependence on the quenching 
factors is depicted. As expected gNa is unconstrained, as shown in the left plot by 
the flat ID marginal posterior (cyan solid line) while for qi one might claim for a slight 
preference for values around 0.08 although it is statistically insignificant. The quenching 
factors for inelastic, Z mediated scattering result to be less constraint than the elastic 
spin- independent case [87]. The right panel illustrates the correlation between {6, rriDM} 
and the quenching factor on Iodine. There is a clear dependence on qi for masses between 
3 and 30 TeV and splittings in the range 50-100 keV: smaller value of the quenching 
factor favours smaller splitting and lighter masses. All the remaining region does not 
show a correlation between the model parameters and qj. The small island at masses of 
few GeV is due to scattering on Sodium and therefore correlated to gNa- 

The dependence on the astrophysical observables and NFW density profile for 
inelastic scattering is shown in figure 13 with the 3D marginal posteriors for {5, mDn} 
and a third parameter direction vq, Wesc and p©. The DAMA signal favours the high tail 
of the velocity distribution from the central panel, where the larger values of fesc dxe 
preferred. From the left and right panel we see that in the 'croissant '-shaped region the 
internal parts are due to circular velocity below vq and DM density close to 0.2 GeVcm"'^. 
The increase of vq and favours instead the outer parts of the region, in example very 
large mass splitting ~ 190 keV and masses around the TeV scale. 

An analogous behavior holds for velocity distributions arising from the cored 
isothermal halo. The main difference is that this latter prefers in particular the very 
high end of the observationally allowed escape velocities. In table 3 the preferred values 



Asymmetric Inert Doublet DM 



27 




50 100 150 200 

6 (keV) 



Figure 12. Inference for DAMA assuming the NFW density profile for the DM 
halo. Left: ID marginal posterior pdf for the quenching factor of Sodium and Iodine, 
as labeled. Right: 3D marginal posterior for {5, mDM,9i}, where the gi direction is 
represented by the colour code. 






Figure 13. Inference for DAMA assuming the NFW density profile for the DM halo. 
Left: 3D marginal posterior for {S,mi)MiVo}i where the vq direction is represented by 
the colour code. Center, right: Same as left, but for {S^muu, Vcsc} and {S^muu, Pq} 
respectively. 



for the astrophysical observables are indicated for both the NFW dark matter profile 
and the cored isothermal one. We underline that the difference in the preferred values 
of fesc will play a role in case of inelastic scattering, even if the statistical significance 
in the difference of the preferred fesc values is small. Indeed in figure 14 we show all 
the experimental constraints and the DAMA region in a single plot, on the left for 
NFW profile and on the right for isothermal cored DM density profile. Firstly we note 
that the NFW profile favours larger splitting for fitting DAMA with respect to the 
isothermal profile and secondly the low mass region is larger. The exclusion limits for 
CDMS (blue dashed), XelOO (pink dot-dashed) and CRESST (black dotted) are 90s% 
confidence intervals and all the region on the left of the curve is excluded. As a general 
remark the NFW prefers smaller Vesc, namely the tail of the velocity distribution is 
constituted by less high speed particles than the isothermal one: there is less room for 
the detectors to be sensitive to inelastic scattering and therefore the exclusion limits 
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Table 3. ID posterior modes and 90% credible intervals for the circular velocity vq, 
escape velocity Veac, and the local DM density for DM density profiles considered 
in this work. 



Vq (km s ^) Vesc (km s ^) Pq (GeV cm ^) 



Cored Isothermal 








DAMA 




629122 


n oi +0.05 
^■'^^-0.03 


CDMS 

XelOO 

CRESST 


211+26 


6291^? 
6281^^ 
628123 


0.31 ±0.04 

si+°-°^ 

0.311^! 


NFW 








DAMA 
CDMS 
XelOO 
CRESST 


220+11 
221+i 
220ttl 


558+20 
558li^ 
557li 
558li7 


0.381°:!^ 

o.38lg:i^ 

0.38lg:it 
0.38lg:i^ 



are less constraining. Regarding the NFW profile, up to masses of 80 GeV the most 
constraining upper bounds is CDMS, then leaving the place to XelOO that excluded all 
DM masses above 316 GeV. The trend for the isothermal profile is the same, except 
that CDMS and XelOO intercepts at 50 GeV and all masses above 251 are excluded by 
XelOO. The transparent region below 50 GeV is excluded by LEPII constraint on the 
Z decay width, while the orange region above 56 TeV is excluded requiring unitarity of 
the S matrix [103]. The CRESST upper bound is comparable to the XelOO one up to 
masses of 300 GeV, because even though it has a much smaller total exposure the W 
is heavier than Xe. For larger masses the effect of the larger total exposure of XelOO 
dominates. 

In terms of SDDM, an asymmetric candidate is completely excluded: the parameter 
space that survives the Xo ~ Xo oscillation bound is severely disfavoured by XelOO, for 
both DM density profiles. In case of the inert doublet model as standard thermal relic 
there is still room up from 45 GeV to 300 GeV for a NFW profile and 250 GeV for an 
isothermal halo. It has been shown that this mass range provides the correct WMAP 
relic abundance thanks to three body annihilation channels [104, 105] and is in the reach 
of LHC [106, 107] . On the contrary, the asymmetric fermionic doublet is a good DM 
candidate up to 200 or 300 GeV depending on the DM density profile, in addition of 
satisfying the DM to baryon ratio. 

6. Conclusions 

We proposed a simple extension of the SM model to explain the observed ratio 
^dm/^b ~ 5 as given by WMAP. We extended the SM by including two heavy triplet 
scalars whose partial decay to SM leptons and inert (odd under a Z2 symmetry) doublet 
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Figure 14. Left: 2D credible regions for the individual experimental bounds and 
regions assuming the NFW DM density profile and marginalizing over the astrophysical 
uncertainties, combined in a single plot. For DAMA (shaded) we show the 90% and 
99% contours. The 90s% contours are given by respectively the pink dot-dash curve 
for XelOO, the dashed blue line for CDMS and the dark green one denotes CRESST. 
The region below = 2 TeV is excluded by xo — Xo oscillation (horizontal red solid 
line), while the orange/dark grey region is excluded by unitarity bound and below the 
dashed gray line by the LEP constraints on the Z decay width. Right: same as left but 
for the isothermal cored DM density profile. 



scalars (x), or vector like fermions (ip), could explain a common origin of asymmetric 
dark matter and visible matter via leptogenesis route. Moreover, the induced vev of 
the triplets also gave rise to neutrino masses, as required by the oscillation experiments, 
via the type-II seesaw mechanism. Thus a triple unification of asymmetric dark matter, 
leptogenesis and neutrino masses could be achieved. 

We studied the relevant annihilation and scattering processes that arise in the 
model. The asymmetry in case of inert scalar (x) doublet dark matter (SDDM) gets 
strongly depleted by the contact annihilation process XX ^ HH mediated via A5 
coupling. Therefore, the survival of the asymmetry in case of inert SDDM required 
A5 < 10~^. Besides that we showed that A5 ~ 10~^ is required for the annual modulation 
signal at DAMA while restoring the asymmetry. On the other hand, the inert fermion 
{if)) doublet dark matter (FDDM) does not under go any further depletion of asymmetry 
in comparison to leptons. A strong constraint arose on the mass scale of inert SDDM 
from the rapid oscillation between xo and its complex conjugate xq. Below EW phase 
transition the fast oscillation between xo and xo depletes the asymmetry strongly. 
Therefore, the survival of the asymmetry in case of SDDM led to its mass > 2 TeV 
so that it freezes out before it begins to oscillate. On the other hand, in case of inert 
FDDM, the survival of asymmetry does not depend on its mass apart from the LEP 
constraint that > Mz/2. Hence a 0(100) GeV dark matter is allowed. 

We then numerically solved the relevant Boltzmann equations to estimate the 
efficiency factors of DM and lepton in either scenarios, for a fixed scalar triplet mass of 
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Ml = 10^° GcV. The model parameter space has been systematically investigated via 
MCMC techniques. We have singled out the preferred regions in the parameter space 
that lead to a successful leptogenesis and to an asymmetric DM, namely satisfying 
fl-Du/^b and fib/ury ratios. We showed that: 

(i) dark matter, irrespective of SDDM or FDDM, masses up to 0{ TeV) can fulfill the 
requirement of Q-Du/^b, 

(ii) for observed BAU and asymmetric dark matter large Bl and small Bdm are 
preferred. In particular for Bl 1 and -Bdm ^ 10~^ the efficiency ratio tjbm/vl 
approaches its maximum value. 

The survival of asymmetry in the dark sector leads to inelastic dark matter because 
the elastic scattering is subdominant in both (SDDM and FDDM) cases. In case of 
SDDM the small coupling A5 ~ 10^^ gave rise to a mass difference between the excited 
state and ground state of DM to be (9(100) keV. On the other hand, in case of inert 
FDDM, the C(IOO) keV mass difference between the ground state and excited state of 
DM is provided by its Majorana mass induced by the triplet scalar. By performing 
a bayesian analysis we found that an asymmetric SDDM of mass larger than 2 TeV 
is strongly disfavoured by the XcnonlOO data while an asymmetric FDDM of mass 
(9(100) GcV is suitable to explain DAMA annual modulation signal while passing the 
latest constraint from XenonlOO experiment. 
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Appendix A. Bayesian inference and MCMC techniques 

In the analysis of the data X the inference of the posterior probability density as a 
function of the parameters, V{9\X), is constructed invoking B ayes' theorem: 

v{e\x)de(x c{x\e) ■ 7i{e)de , (a.i) 

where C{X\9) is the likelihood function and 7i{6) denotes the probability density on the 
parameter space 9 prior to observing the data X. The posterior pdf represents our state 
of knowledge about the parameters after taking into account the information contained 
in the data, and has an intuitive and straightforward interpretation in that Jy V{9\X)d9 
is the probability that the true value of 9 lies in the volume V. 

Since the prior pdf is independent of the data, it needs to be chosen according 
to one's belief and is thus inherently subjective. In the often encountered situation in 
which no unique theoretically motivated prior pdf can be derived, one may wish to use 
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one which does not favour any parameter region in particular. A common choice is the 
uniform or top-hat prior 

" 1 0, otherwise, ^^"^^ 

if the general order of magnitude of the parameter is known. Here, the limits O^i^ and 
^max should be chosen such that they are well beyond the parameter region of interest. 
If even the order of magnitude is unknown, one may want to choose a uniform prior in 
log 9 space instead, 

n a f dlog6', if 6'min < 6* < 6'max, /.ON 

.,„,(log^)dlog^=| ^^^^^^.^^^ (A.3) 

which is equivalent to a d6/6 prior in 6 space. Note that because the volume element 
d^ is in general not invariant under a parameter transformation f : 6 ^ 6' , a uniform 
prior pdf on 6 does not yield the same probabilities as a uniform prior pdf on 6' unless 
the mapping / is linear. The same is also true for the posterior probabilities, i.e., 
V{e\X)de ^ V{e'\X)de' m general. 

While the posterior pdf technically contains all the necessary information for the 
interpretation of the data, the fact that it is a function in the A'"-dimensional space 
of parameters makes it difficult to visualise if iV > 2. Being a probability density, 
its dimensionality can be easily reduced by integrating out less interesting {nuisance) 
parameter directions ipi, yielding an n-dimensional marginal posterior pdf, 

(X J di;i...d7pmV{9i,...,9n,ipi...,ipm\X), (A.4) 

which is more amenable to visual presentation if n = 1,2, and can be used to construct 
constraints on the remaining parameters. 

We employ a modified version of the generic Metropolis-Hastings sampler [108, 109] 
included in the public MCMC code CosmoMC [110, 111], to sample the posterior over the 
full parameter space. Each point Xi+i obtained from a Gaussian random distribution 
(the so called proposal density) around the previous point is accepted to be the next 
element of the chain with probability: 

The resulting chains are analysed with an adapted version of the accompanying 
package GetDist, supplemented with matlab scripts from the package SuperBayeS [112, 
113]. One- or two-dimensional marginal posterior pdfs are obtained from the chains 
by dividing the relevant parameter subspace into bins and counting the number of 
samples per bin. An x% credible interval or region containing x% of the total volume 
of ■Pmar IS then constructed by demanding that Pmar at any point inside the region be 
larger than at any point outside. In the one-dimensional credible interval thus 

constructed corresponds to the Minimal Credible Interval of [114]. Provided the data are 
sufficiently constraining the marginal posterior typically exhibits very little dependence 
on the choice of prior. For data that can only provide an upper or a lower bound 



P(x.^O = min(l, . (A.5) 
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on a parameter (or no bound at all) however, the properties of the inferred posterior 
and the boundaries of credible regions can vary significantly with the choice of prior as 
well as its limits ^^min and ^^max, making an objective interpretation of the results rather 
difficult. This is in particular the case for the inference of credible regions subject to 
direct detection data. 

Indeed in the case of XelOO, CDMS and CRESST we construct intervals on the 
volume of the marginal posterior in 5'-space Vraar{S\X), where S is the expected WIMP 
signal, using a uniform prior on S with a lower boundary at zero. An x% upper 
bound thus constructed has a well-defined Bayesian interpretation that the probability of 
S < S^is x%. The limit is then mapped onto the {moM, (5}-plane by identifying those 
combinations of ttidm and S with 'Pmar('^DM, ^|-^) = VmariSxlX). An x% contour, which 
will be denoted Xs%, computed in this manner has the property of being independent 
of the choice of prior boundaries for ttidm and S. Its drawback, however, is that it has 
no well-defined probabilistic interpretation in {uidm, S}-spa.ce. 

Appendix B. Standeird Model Halo results for ineleistic scattering 

For sake of completeness we show in figure Bl the preferred region for DAMA modulation 
assuming the standard maxwellian halo with fixed astrophysical observablcs as DM 
velocity distribution, Eq. 5.2. The astrophysical variables are fixed at their preferred 
values, sec section 5.1. In the same plot we add the exclusion limits of CDMS, XelOO and 
CRESST. The trend is analogous to figure 14, even though the DAMA region is smaller, 
since there a no volume effect due to marginalization over the astrophysical uncertainties. 
The region is also smaller because of Vesc — 544 kms~^, which limits the contribution of 
particles in the very high tail of the distribution. At the same time the exclusion bounds 
are tighter than in NEW and isothermal case, again because the reduce parameter space. 
The results regarding the viability of the DM candidate are practically unchanged with 
respect to the case with marginalization over the astrophysics: an asymmetric scalar 
DM with ^ 2TeV as explanation of the DAMA signal is completely disfavoured by 
XelOO, while a fermionic component is compatible at 90% C.L. up to a mass of 200 
GeV. 

References 

[1] Komatsu E et al. (WMAP Collaboration) 2011 Astrophys.J.Suppl. 192 18 {Preprint 1001.4538) 
[2] Kohri K, Mazumdar A and Sahu N 2009 Phys.Rev. D80 103504 {Preprint 0905.1625) 
[3] Nussinov S 1985 Phys.Lett. B165 55 

[4] Barr S M, Chivukula R and Farhi E 1990 Phys.Lett. B241 387-391 

[5] Dodelson S, Greene B R and Widrow L M 1992 Nucl.Phys. B372 467-493 

[6] Kaplan D B 1992 Phys.Rev. Lett. 68 741-743 

[7] Kuzmin V A 1998 Phys.Part.Nucl. 29 257 265 {Preprint hep-ph/9701269) 

[8] Fnjii M and Yanagida T 2002 Phys.Lett. B542 80 88 {Preprint hep-ph/0206066) 

[9] Oaknin D H and Zhitnitsky A 2005 Phys.Rev. D71 023519 {Preprint hep-ph/0309086) 



Asymmetric Inert Doublet DM 



33 




Figure Bl. 2D credible regions for tlie individual experimental bounds and regions 
assuming the SMH, combined in a single plot. For DAMA (shaded) we show the 
90% and 99% contours. The 90s% contours are given by respectively the pink dot- 
dash curve for XelOO, the dashed blue line for CDMS and the dark green one denotes 
CRESST. The region below AI^ — 2 TeV is excluded by xo — Xo oscillation (horizontal 
red solid line), while the orange/dark grey region is excluded by unitarity bound 
and below the dashed gray line by the LEP constraints on the Z decay width. The 
astrophysical observables are fixed at vq, Vcsc and Pq. 



[10] Hooper D, March-Russell J and West S M 2005 Phys.Lett. B605 228-236 {Preprint 
hep-pli/0410114) 

[11] Kitano R and Low I 2005 Phys.Rev. D71 023510 {Preprint hep-ph/0411133) 
[12] Cosme N, Lopez Honorez L and Tytgat M H 2005 Phys.Rev. D72 043505 {Preprint 
hep-ph/0506320) 

[13] Farrar G R and Zaharijas G 2006 Phys.Rev. Lett. 96 041302 {Preprint liep-pli/05 10079) 

[14] Roszkowski L and Seto O 2007 Phys.Rev. Lett. 98 161304 {Preprint hep-ph/0608013) 

[15] Kaplan D E, Luty M A and Zurek K M 2009 Phys.Rev. D79 115016 {Preprint 0901.4117) 

[16] Kohri K, Mazumdar A, Sahu N and Stephens P 2009 Phys.Rev. D80 061302 {Preprint 0907 . 0622) 

[17] An H, Chen S L, Mohapatra R N and Zhang Y 2010 JHEP 1003 124 {Preprint 0911.4463) 

[18] Frandsen M T and Sarkar S 2010 Phys.Rev. Lett. 105 011301 {Preprint 1003.4505) 

[19] Feldstein B and Fitzpatrick A 2010 JCAP 1009 005 {Preprint 1003.5662) 

[20] An H, Chen S L, Mohapatra R N, Nussinov S and Zhang Y 2010 Phys.Rev. D82 023533 {Preprint 
1004.3296) 

[21] Cohen T, Phalen D J, Pierce A and Zurek K M 2010 Phys.Rev. D82 056001 {Preprint 1005 . 1655) 

[22] Shelton J and Zurek K M 2010 Phys.Rev. D82 123512 {Preprint 1008.1997) 

[23] Davoudiasl H, Morrissey D E, Sigurdson K and Tulin S 2010 Phys.Rev. Lett. 105 211304 {Preprint 



Asymmetric Inert Doublet DM 



34 



1008.2399) 

[24] Gu P H, Lindner M, Sarkar U and Zhang X 2010 {Preprint 1009 . 2690) 

[25] Blennow M, Dasgupta B, Fernandez-Martinez E and Rius N 2011 JHEP 1103 014 {Preprint 
1009.3159) 

[26] McDonald J 2010 {Preprint 1009 . 3227) 

[27] Hall L J, March-Russell J and West S M 2010 {Preprint 1010.0245) 

[28] Dutta B and Kumar J 2011 Phys.Lett. B699 364-367 {Preprint 1012.1341) 

[29] Haba N, Matsumoto S and Sato R 2011 {Preprint 1101.5679) 

[30] Falkowski A, Rudcrman J T and Volansky T 2011 JHEP 1105 106 {Preprint 1101.4936) 

[31] Chun E J 2011 JHEP 1103 098 {Preprint 1102.3455) 

[32] Kang Z, Li J, Li T, Liu T and Yang J 2011 {Preprint 1102.5644) 

[33] Heckman J J and Rey S J 2011 JHEP 06 {Preprint 1102.5346) 

[34] Frandsen M T, Sarkar S and Schmidt-Hoberg K 2011 {Preprint 1103.4350) 

[35] Buckley M R 2011 {Preprint 1104.1429) 

[36] Iminniyaz H, Drees M and Chen X 2011 {Preprint 1104.5548) 

[37] March-Russell J and McCuUough M 2011 {Preprint 1106.4319) 

[38] Del Nobile E, Kouvaris C and Sannino F 2011 Phys.Rev. D84 027301 {Preprint 1105.5431) 

[39] McDermott S D, Yu H B and Zurek K M 2011 {Preprint 1103.5472) 

[40] Kouvaris C and Tinyakov P 2011 {Preprint 1104.0382) 

[41] Deshpande N G and Ma E 1978 Phys. Rev. D18 2574 

[42] Ma E 2006 Phys.Rev. D73 077301 {Preprint hep-ph/0601225) 

[43] Barbieri R, Hall L J and Rychkov V S 2006 Phys.Rev. D74 015007 {Preprint hep-ph/0603188) 

[44] Lopez Honorez L, Nezri E, Oliver J F and Tytgat M H 2007 JCAP 0702 028 {Preprint 
hep-ph/06 12275) 

[45] Hambye T and Tytgat M H G 2008 Phys. Lett. B659 651-655 {Preprint 0707.0633) 

[46] Magg M and Wetterich C 1980 Phys.Lett. B94 61 

[47] Cheng T and Li L F 1980 Phys.Rev. D22 2860 

[48] Schechter J and Valle J 1980 Phys.Rev. D22 2227 

[49] Gelmini G and Roncadelli M 1981 Phys.Lett. B99 411 

[50] Lazarides G, Shafi Q and Wetterich C 1981 Nucl.Phys. B181 287 

[51] Mohapatra R N and Senjanovic G 1981 Phys.Rev. D23 165 

[52] Fukuda S et al. (Super-Kamiokande Collaboration) 2001 Phys. Rev. Lett. 86 5656-5660 {Preprint 
hep-ex/0103033) 

[53] Ahmad Q et al. (SNO Collaboration) 2002 Phys. Rev. Lett. 89 011302 {Preprint nucl-ex/0204009) 

[54] Eguchi K et al. (KamLAND CoUaboration) 2004 Phys. Rev. Lett. 92 071301 {Preprint 
hep-ex/0310047) 

[55] Nezri E, Tytgat M H G and Vcrtongcn G 2009 JCAP 0904 014 {Preprint 0901.2556) 

[56] Tucker-Smith D and Weiner N 2001 Phys. Rev. D64 043502 {Preprint hep-ph/0101138) 

[57] Bernabei R, Belli P, Cappella F, CeruUi R, Dai C et al. 2010 Eur.Phys.J. C67 39-49 {Preprint 
1002.1028) 

[58] March-Russell J and West S M 2004 Phys. Lett. B593 181 188 {Preprint liep-ph/0403067) 

[59] March-Russell J, McCabe C and McCullough M 2008 {Preprint 0812.1931) 

[60] Arrenberg S, Baudis L, Kong K, Matchev K T and Yoo J 2008 Phys. Rev. D78 056002 {Preprint 
0805.4210) 

[61] Cui Y, Morrissey D E, Poland D and Randall L 2009 {Preprint 0901.0557) 

[62] Alves D S M, Behbahani S R, Schuster P and Wacker J G 2009 {Preprint 0903.3945) 

[63] Finkbeiner D P, Slatyer T R, Weiner N and Yavin I 2009 {Preprint 0903. 1037) 

[64] Finkbeiner D P, Lin T and Weiner N 2009 {Preprint 0906 . 0002) 

[65] Farina M, Kadastik M, Pappadopulo D, Rata J, Raidal M et al. 2011 {Preprint 1104.3572) 

[66] Schwetz T and Zupan J 2011 * Temporary entry * {Preprint 1106.6241) 

[67] Arina C, Ling F S and Tytgat M H 2009 JCAP 0910 018 {Preprint 0907.0430) 



Asymmetric Inert Doublet DM 



35 



[68] Ahmed Z et al. (CDMS Collaboration) 2010 [Preprint 1012.5078) 
[69] Angloher G et al. 2008 [Preprint 0809 . 1829) 

[70] Aprile E et al. (XENONIOO Collaboration) 2011 [Preprint 1104.3121) 

[71] Ma E and Sarkar U 1998 Phys. Rev. Lett. 80 5716-5719 [Preprint hep-ph/9802445) 

[72] Hambyc T, Ma E and Sarkar U 2001 Nucl.Phys. B602 23-38 [Preprint hep-ph/0011192) 

[73] Sahu N and Sarkar U 2007 Phys.Rev. D76 045014 [Preprint hep-ph/0701062) 

[74] McDonald J, Sahu N and Sarkar U 2008 JCAP 0804 037 [Preprint 0711.4820) 

[75] Majcc S K and Sahu N 2010 Phys.Rev. D82 053007 [Preprint 1004.0841) 

[76] Lundstrom E, Gustafsson M and Edsjo J 2008 [Preprint 0810.3924) 

[77] Hambye T and Senjanovic G 2004 Phys. Lett. B582 73-81 [Preprint hep-ph/0307237) 

[78] Sahu N and Sarkar U 2006 Phys.Rev. D74 093002 [Preprint hep-ph/0605007) 

[79] Harvey J A and Turner M S 1990 Phys.Rev. D42 3344 3349 

[80] Hambye T, Raidal M and Strumia A 2006 Phys.Lett. B632 667 674 [Preprint hep-ph/0510008) 

[81] Chun E J and Scopel S 2007 Phys.Rev. D75 023508 [Preprint hep-ph/0609259) 

[82] Clesse S, Ringeval C and Rocher J 2009 Phys.Rev. D80 123534 [Preprint 0909.0402) 

[83] de Austri R R, Trotta R and Roszkowski L 2006 JHEP 0605 002 [Preprint hep-ph/0602028) 

[84] Loredo T J 1995 81-142 Available from http://bayes.wustl.edu/gregory/articles.pdf 

[85] Trotta R 2008 Contemp.Phys. 49 71-104 [Preprint 0803.4089) 

[86] HINKLEY D V 1969 Biometrika 56 635-639 

[87] Arina C, Hamann J and Wong Y Y 2011 [Preprint 1105.5121) 

[88] Helm R H 1956 Phys. Rev. 104 1466-1475 

[89] Lewin J D and Smith P F 1996 Astropart. Phys. 6 87-112 

[90] Navarro J F, Frenk C S and White S D M 1997 Astrophys. J. 490 493-508 [Preprint 
astro-ph/9611107) 

[91] Binney J and Tremaine S 1987 Galactic Dynamics (Princeton University Press) 
[92] Aalseth C, Barbeau P, Colarcsi J, Collar J, Leon J et al. 2011 [Preprint 1106.0650) 
[93] Lebedenko V N e< al. 2008 [Preprint 0812.1150) 

[94] Bernabei R et al. (DAMA) 2008 Eur. Phys. J. C56 333-355 [Preprint 0804.2741) 
[95] Bernabei R et al. 1996 Phys. Lett. B389 757-766 

[96] Chagani H. Majcwski P, Daw E, Kudryavtsev V and Spooner N 2008 JINST 3 P06003 [Preprint 
0806.1916) 

[97] Reid M, Menten K, Zheng X, Brunthaler A, Moscadelli L et al. 2009 Astrophys. J. 700 137-148 
[Preprint 0902.3913) 

[98] Gillcsscn S, Eisenhauer F, Trippe S, Alexander T, Genzel R et al. 2009 Astrophys.J. 692 1075- 

1109 [Preprint 0810.4674) 
[99] Smith M C, Ruchti G, Helmi A, Wyse R, Fulbright J et al. 2007 Mon.Not.Roy.Astron.Soc. 379 

755-772 [Preprint astro-ph/0611671) 
[100] Dehnen W and Binney J 1998 Mon. Not. Roy. Astron. Soc. 298 387-394 [Preprint 

astro-ph/9710077) 

[101] Weber M and de Boer W 2010 Astron.Astrophys. 509 A25 [Preprint 0910.4272) 
[102] Salucci P, Nesti F, Gentile G and Martins C 2010 Astron.Astrophys. 523 A83 [Preprint 
1003.3101) 

[103] Griest K and Kamionkowski M 1990 Phys. Rev. Lett. 64 615 

[104] Lopez Honorez L and Yaguna C E 2010 .JHEP 1009 046 [Preprint 1003.3125) 

[105] Lopez Honorez L and Yaguna C E 2011 JCAP 1101 002 [Preprint 1011 . 1411) 

[106] Dolle E, Miao X, Su S and Thomas B 2010 Phys.Rev. D81 035003 [Preprint 0909.3094) 

[107] Miao X, Su S and Thomas B 2010 Phys.Rev. D82 035009 [Preprint 1005.0090) 

[108] Metropolis N, Roscnbluth A W, Rosenbluth M N, Teller A H and Teller E 1953 Journal of 

Chemical Physics 21 1087-1092 
[109] Hastings W 1970 Biometrika 97-109 

[110] Lewis A and Bridle S 2002 Phys. Rev. D66 103511 [Preprint astro-ph/0205436) 



Asymmetric Inert Doublet DM 



36 



[111] Lewis A and Bridle S CosmoMC Notes http://cosmologist.info/notes/CosmoMC.pdf 
[112] de Astrui R R, Trotta R and Feroz F SuperBayeS package http://www.superbayes.org/ 
[113] Trotta R, Feroz F, Hobson M P, Roszkowski L and Ruiz de Austri R 2008 JHEP 0812 024 
{Preprint 0809 . 3792) 

[114] Hamann J, Hannestad S, Raffelt G and Wong Y Y 2007 JCAP 0708 021 {Preprint 0705.0440) 



